Conditional generation of arbitrary multimode entangled states of light 

with linear optics 
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We propose a universal scheme for the probabilistic generation of an arbitrary multimode entan- 
gled state of light with finite expansion in Fock basis. The suggested setup involves passive linear 
optics, single photon sources, strong coherent laser beams, and photodetectors with single-photon 
resolution. The efficiency of this setup may be greatly enhanced if, in addition, a quantum memory 
is available. 
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I. INTRODUCTION 

The generation of nonclassical states of light is one of 
the primary research areas in quantum optics. In partic- 
ular, the preparation of entangled states of light has at- 
tracted a considerable amount of attention recently, since 
these states have been identified as a key resource for 
quantum information processing [lj. Other interesting 
applications of the entangled states include,for instance, 
ultra- high precision measurements EH 01 13 an d quan- 
tum optical lithography 0, Q • 

The range of interactions between light fields that are 
experimentally accessible and feasible is rather limited, 
which restricts the class of quantum states of the op- 
tical field that can be generated in the lab. However, 
this class can be significantly extended if one consid- 
ers probabilistic generation schemes, whose success is 
conditioned on the detection of a particular outcome 
of a measurement performed on some ancilla system. 
Schemes for probabilistic preparation of Fock states , 
arbitrary superpositions of Fock states of single-mode 
field Gl[ll| , and superpositions of classically distinguish- 
able states (Schrodinger-cat-like states) [13, 13], have 
been found. Several setups for the generation of two- 
mode TV-photon path-entangled states have been sug- 
gested [H E El |H 13. Recently, the experimen- 
tal conditional preparation of a single-photon Fock state 
with negative Wigner function has been reported [l9| . 

However, all the above mentioned schemes are still 
quite restrictive since they are capable to prepare only 
single-mode states or only a particular class of two-mode 
states (such as the iV-photon states). In this paper, we 
design a truly universal scheme that can be used to prob- 
abilistically generate an arbitrary multimode entangled 
state of light, provided that each mode does not contain 
more than N photons, where N is an arbitrary but finite 
integer. The resources required for the present scheme 
comprise the passive linear optical elements (beam split- 
ters and phase shifters), single photon sources, strong 
coherent laser beams, and photodetectors. It was shown 
recently that these resources are sufficient for perform- 
ing universal quantum computing pfl Eil I2II23T I2I |25| . 



Our proposal provides another illustration of the surpris- 
ing versatility and power of such an approach, relying 
only on linear optics and single photons. In this con- 
text it is worth mentioning that, very recently, Clausen 
et al. suggested a scheme for approximate conditional 
implementation of general single- or multi-mode opera- 
tors acting on states of traveling optical field [26j. The 
scheme described in [2(| is also based on linear optics 
and may be in principle used for (approximate) quantum 
state preparation. We emphasize that the approach pro- 
posed in the present paper conceptually differs from that 
of Ref. [26|, We are primarily interested in exact state 
preparation and our setup is specifically tailored for that 
purpose. 

We will first explain all the essential features of the 
state-preparation procedure on the simplest yet nontriv- 
ial example of two-mode entangled state that is formed 
by a superposition of two product states, 
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where Af is a normalization factor. Then we will gener- 
alize the preparation procedure and design a scheme for 
the generation of arbitrary multimode entangled states of 
light. We note in passing that even the simple entangled 
states Q may find applications in the investigations of 
the foundations of quantum mechanics. For instance, it 
was shown recently that the states of the form Q yield 
a strong violation of Bell-CHSH inequalities when Alice 
and Bob perform balanced homo dyne measurements fol- 
lowed by an appropriate binning |27| . 

The paper is structured as follows. In Section II, we 
propose a method to probabilistically generate the two- 
mode entangled state with the help of passive linear 
optics, single photon sources, strong coherent laser fields 
and photodetectors. In section III we shall extend our 
scheme to allow for more than two modes and more than 
two terms in the superposition. We will thus design a 
universal device that enables one to conditionally pre- 
pare an arbitrary multimode state of light. Finally, the 
conclusions are drawn in Section IV. 
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II. CONDITIONAL GENERATION OF 
TWO-MODE ENTANGLED STATES 

Let us first note that, in order to simplify the notation, 
we shall omit the normalization factors in front of the 
quantum states in what follows. 

Since the state \tp}AB is rather complicated, we divide 
its preparation into several steps. The main simplifica- 
tion stems from the observation that the entangled state 
itTll can be prepared by means of entanglement swapping 
|28l I2H l30| if we possess the two three-mode states 

\4>)a = \f)A 1 \V)A 2 + \g) Al \H) A2 , (2) 
W)b = \f')B 1 \V) B2 + \g')B 1 \H)B 2 , (3) 

where \V) and \H) denote the state of a single photon 
in spatial mode A2 (or B2) that is polarized vertically or 
horizontally, respectively. As shown in Fig. 1, we easily 
obtain the state from the state I^aI^') b if we project 
the single photons in spatial modes Ai and B2 onto the 
maximally entangled Bell state 
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The generation of state Q thus boils down to the prepa- 
ration of the entangled state J2J)- [State can be pre- 
pared similarly]. This latter task will be accomplished 
with the help of a quantum non-demolition (QND) mea- 
surement of a single photon, which creates entanglement. 



A. QND measurement 

The proposed setup is depicted in Fig. 2. The input 
polarization modes IV and 1H are prepared in the pure 
single-mode states |/) and \g) that contain no more than 
N + 1 photons, 
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where \n) denotes the n-photon Fock state. The states 
© are, of course, closely related to the states |/) = 
E^=o/«l n ) and Iff) = T,n=o9n\n) appearing in Eq. ©. 




FIG. 1: Schematic of the setup for the generation of the en- 
tangled state Q via entanglement swapping. The two modes 
emerging from the QND measurement devices are combined 
on a balanced beam splitter and the swapping succeeds when 
each photodetector PD detects exactly one photon. The boxes 
denoted Alice and Bob correspond to the setup depicted in 
Fig. 2 and WP denotes a waveplate. 
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FIG. 2: The setup for the conditional generation of the entan- 
gled state The polarizing beam splitter PBS\ transmits 
vertically polarized photons and reflects horizontally polar- 
ized ones. The box labeled as QND performs the quantum 
non-demolition measurement of the photon number in the 
spatial mode while preserving the polarization state of the 
photon, see Fig. 3. The polarizing beam splitter PBS2 par- 
tially reflects and partially transmits both vertically and hor- 
izontally polarized photons. BSi is a usual beam splitter and 
PD denotes photodetector with single-photon sensitivity. 



The exact relationship between them will be specified 
later. Note that it was shown by Dakna et al. that 
the single-mode states (0 can be probabilistically gener- 
ated with the help of passive linear optics, single photon 
sources, strong coherent laser pulses and photodetectors 
with single photon sensitivity. 

As shown in Fig. 2, the input beam containing the 
states I/) and \g) in the V and H polarization modes 
impinges on a beam splitter BS\ with transmittance t\ 
and reflectance n. The purpose of this beam splitter is 
to separate a single photon from one of the input states 
and thus create entanglement. Of course, BS\ may split 
more than one photon or no photon at all. Therefore we 
must verify that there is exactly a single photon present 
in the output spatial mode 2, without disturbing its po- 
larization state. Such QND measurements of a single 
photon have been recently thoroughly discussed in Ref. 
[31| where several schemes have been proposed and ana- 
lyzed. 

One possible method relies on teleportation HH . 
With the help of the conditional C-NOT gate for pho- 
tonic qubits j2£| we have to first prepare the singlet state 
\VH)-a — \HV)zi in the auxiliary modes 34. Then we 
teleport the polarization state of a photon in mode 2 onto 
the photon in mode 4 by performing a Bell measurement 
on modes 2 and 3. This measurement should be carried 
out with two detectors that are capable to resolve the 
number of photons. If exactly two photons are detected 
in coincidence at the two detectors, then we know that 
there was a single photon in mode 2 and its polarization 
state has been coherently transfered onto the state of a 
single photon in mode 4. 
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FIG. 3: Interferometric scheme for a QND measurement of a 
single photon pjlll . BS are balanced beam splitters. PBS de- 
notes polarizing beam splitters that transmit horizontally po- 
larized photons and reflect vertically polarized photons. Rgo 
is a waveplate that rotates polarization state by 90 degrees. 



Kok et al. also proposed QND measurement schemes 
that do not require a-priori entanglement [3]]]. For our 
purposes it suffices to employ the simple interferomet- 
ric scheme depicted in Fig. 3 which performs a partial 
QND photon number measurement. The coincident de- 
tection of a single photon in detectors PD\ and PD 2 
indicates that there was at least one photon in the in- 
put mode. Moreover, if there was exactly a single photon 
at the input, then its polarization state is unperturbed 
by the measurement. The probability of successful QND 
measurement of a single photon provided that there is a 
single photon in the spatial mode 2 is equal to 1/8. In 
the preparation scheme shown in Fig. 1 the two QND 
measurements are followed by the entanglement swap- 
ping that succeeds only if each photodetector in Fig. 1 
detects exactly one photon. In this way we select the 
events when there is exactly a single photon in each mode 
A 2 and B 2 . 

The state after the QND measurement can be written 
in the form 

I^qnd) = Bi|/)iv B \g)iH \V)t+Bo\f)iv \B) a , 

(6) 

where the non-unitary operators Bq and B\ describe re- 
moval of none or a single photon at BSi, respectively, 
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B 1 = ^V^+T^n|n)(n + l|. 



(7) 
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B. Quantum erasure 

The QND measurement has thus created the entangle- 
ment. Still, Bj\f) and Bj\g) are states of two different 



polarization modes while |/) and \g) in Eq. are states 
of a single mode. We achieve this by erasing the infor- 
mation about the polarization. This is done with the 
help of a polarizing beam splitter PBS2, see Fig. 2. The 
beam splitter is rotated such that it combines the V and 
H polarizations. If we associate the annihilation opera- 
tors aiv and am with the modes IV and 1H, then the 
transformation carried out by PBS2 is given by 



aix = r 2 aiv + ham, 
air = r 2 a 1H - t 2 di V , 



(8) 



where the polarization modes IX and 1Y are spatially 
separated at the output of PBS 2 as indicated in Fig. 
2. The erasing succeeds if and only if the detector PD 
placed on the 1Y mode does not detect any photon, which 
ensures that all photons have been transmitted to the lin- 
early polarized output mode IX. 

Since the input states do not contain more than N+l 
photons each and since a single photon was subtracted 
on BSi, the state after erasure of the polarization infor- 
mation reads 



2N+1 



I^pbs) = J2 ifn\n)\V)+g n \n)\H)). 



(9) 



where the output complex amplitudes f n and g n are given 
by the following formula, 
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(10) 



It follows that we can fix the values of output coeffi- 
cients /„ by properly choosing the complex amplitudes 
of the input states. We have the recurrence equations 
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*?rir§<2 k fk+ig n -k, 



(11) 



and a similar formula holds for g n +i- The amplitudes 
fn+i and g n +i are calculated as follows. First, some 
nonzero values for /q and go are chosen. Then one uses 
repeatedly Eq. I jllj l to calculate f n +\ and g n +i from //. 
and <jk, < k < n, and the iterations stop at n = N. 

We have thus almost obtained the state (J2J. How- 
ever, while the effective states |/} and \g) in |JSJ have the 
correct structure in the subspace spanned by the first 
N + l Fock states |0), |1), . . . , \N), they also contain an 
unwanted tail in the subspace spanned by the Fock states 
|JV + 1),...,|2JV + 1). 
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FIG. 4: Quantum scissors. 



should hold, where 
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(14) 



is the probability of the successful quantum-state trun- 
cation that projects onto the subspace spanned by 
|0),|l),...,|iV). 



III. GENERAL SCHEME 



C. Quantum state truncation 

The last step in our scheme is to get rid of that tail. 
We do so with the help of the quantum scissors [3, HE 
l3ri l37l l38l | that will project the state in the mode 1 onto 
the subspace spanned by the first N + 1 Fock states. 
This action is formally described by the projector Hqs = 

^2n=o \ n )( n \- After this truncation, we finally obtain the 
target entangled state © in the output modes 1 and 2 
of our device. 

The quantum scissors have been discussed in detail in 
several recent papers [2J, t&J l2a biZl till ■ The schemes 
proposed there can be implemented with the use of linear 
optics, single photon sources, and photodetectors. The 
setup that we shall consider in the present paper is de- 
picted in Fig. 4. The modes 2 and 3 are prepared in a 
pure iV-photon two-mode entangled state 
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Y,c k \k) 2 \N - k)< 



(12) 



k=0 



It was shown recently that these states may be condi- 
tionally generated with the help of linear optics, single 
photon sources and photodetectors [lj, [Tj, [T|| [13 • We 
then combine the modes 1 and 2 on a balanced beam 
splitter and measure the number of photons in output 
modes 1 and 2 by means of two photodetectors PD\ and 
PD2 with single-photon resolution. The truncation is 
a conditional operation that succeeds when certain par- 
ticular measurement outcome has been detected, say, no 
photons at PD\ and N photons at PD2, which we denote 
as (0,N). Note, that one could also consider other de- 
tection events, but the total number of photons detected 
by the photodetectors should be equal to TV. The quan- 
tum scissors is essentially a quantum state teleportation 
in the subspace of the first N + 1 Fock states, where the 
state 1)1 21) serves as the quantum channel, and the Bell- 
type measurement is carried out with the use of the beam 
splitter in Fig. 4 that couples the modes 1 and 2. 

For our particular choice of the detection event (0, N), 
and assuming a balanced beam splitter, we find that 



QS" 



■y/k\{N-k)\ 



(13) 



So far we have considered the generation of a two- 
mode state that consists of a superposition of two product 
states. In this section we shall generalize our scheme in 
several ways. 



A. M-mode entangled state 

First, we will discuss the preparation of a M-modc 
entangled state which has a form limited to the superpo- 
sition of two product states, namely 

\4') = \fi)\f2)...\fM)+e* e \ 9l }\g 2 }...\gM). (15) 

To prepare this state, we first generate M states 

\<t>j) = \fj)\V) + \gj)\H) (16) 

as before, and then project the M single photons onto 
the maximally entangled state 



\GHZ) = \VV ...V) 



-i0 
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(17) 



This projection can be performed with the Greenberger- 
Horne-Zcilinger (GHZ) state analyzer described by Pan 
and Zeilinger |3Jj |. This analyzer consists of polarizing 
beam splitters and half-wave plates and allows one to 
distinguish two GHZ states of M qubits. 



B. Arbitrary two- mode entangled states 

Our second generalization extends the preparation pro- 
cedure to two-mode states that are superpositions of d 
product states, 



)AB=Y, el63 \fi)^j)B. 
3=1 



(18) 



Obviously, a natural strategy is to prepare this state via 
c?-dimcnsional entanglement swapping from two states of 
the form 



3=1 



)BD = £ \9j)B\l)jD, 
3=1 

(19) 
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FIG. 5: Scheme for conditional generation of the entangled 
state H19fl . The erasing of the spatial information is achieved 
conditionally on no photon being detected at the d — 1 pho- 
todetectors PD. 



anced linear superposition of all d modes 
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(21) 



The complex amplitudes of the states \ fj) — ^ n /j, n |^) 
can then be expressed in terms of the complex amplitudes 
of the input states fj in as follows 



d n/2 



V n 3 + 1 



fj,n } +l TT fk,n k 



n 



'n k \ 



(22) 



where the prime indicates summation over all n — 
(ni, . . . , n d ) satisfying the constraint X)fc=i n k = n. Note 
that the formula (|22|l may be easily inverted and we can 
determine the complex amplitudes fj t k of the input states 
for any given prescribed output states \fj). 



Universal scheme 



where \l)j denotes a state of a single photon in the jth 
auxiliary spatial mode. To prepare the states (|19fl . we 
proceed essentially as in Sec. II. The extended scheme 
is depicted in Fig_. 5. Initially, d spatial modes are pre- 
pared in states \fj). The array of d beam splitters BSj 
extracts, with a certain probability, a single photon. We 
must verify the presence of exactly one single photon in 
total in the modes 1', . . . , dl while preserving the coher- 
ence among these modes by performing a QND measure- 
ment. This can be accomplished by d-dimensional tele- 
portation similarly as discussed in Section II. However, 
this time we need to teleport a qudit encoded as a state 
of a single photon in d spatial modes. Fortunately, it 
was shown recently that the probabilistic gates for quan- 
tum computing with linear optics [23, El H3, El EU- 
that were introduced for qubits, can be easily extended 
to qudits ^3]- This means that all the necessary manip- 
ulations, such as the generation of maximally entangled 
state of two qudits and the Bell measurement, can be 
probabilistically carried out with the resources that we 
assume here. 

After the QND measurement, we must erase the spatial 
information contained in the resulting state 

d / d \ 

j=i \fe=i / 

where Bj are given by Eq. Q. This is done by the 
second array of d — 1 beam splitters indicated by the 
shaded zone in Fig. 5. The preparation is finished by the 
application of the quantum scissors, that truncate the 
Fock state expansion of the output state at N photons. 

Suppose that the output mode after erasing is a bal- 



Finally, we point out that those two generalizations 
can be combined and we can thus generate an arbitrary 
multimode entangled state of light where each mode con- 
tains no more than N photons. Any M-mode state of this 
kind can be written as a superposition of no more than 
d = (N + 1) M product states, 

d 

\^) = Y. eiei \fn)i\h2)2...\fm)M. (23) 

Note that the norm of each product state in the super- 
position 12: il) may be arbitrary. We can conditionally 
generate the state H23H if we first prepare M entangled 
states (k = 1, .. . ,M) 

d 

\4>k) =Y,\fik)\l) jk , (24) 

3=1 

where \l)jk denotes a single photon state of the (j,k)- 
th auxiliary spatial mode. As discussed above, the states 
can be prepared with the help of the scheme depicted 
in Fig. 5. Then, we carry out a joint measurement on all 
auxiliary photons. The state is obtained if we project 
the auxiliary modes onto an M-photon GHZ entangled 
state: 

d 

\GHZ Md ) =^e-^|%|l) i2 ...|l) jM . (25) 

3=1 

If we denote |$Md) = |<^i)|<fe) • ■ ■ \4>m) then we can write 

|V) = (GHZ Md \$Md). (26) 

With the help of the techniques developed in the frame- 
work of the quantum computation with linear optics, the 
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entangled M-photon GHZ state l|25|l can be mapped onto 
a product state of M photons which can then be detected 
by observing an appropriate M-photon coincidence. The 
factorization consists of a sequence of the C-SHIFT gates 
for two qudits. The C-SHIFT gate is defined as follows, 

|i)ifc|i)i'*' - \Vjk\i)r-i, k ', (27) 

where f — j should be calculated modulo d. As shown 
in Ref. 0] this transformation can be (probabilistically) 
performed with the use of single photon sources, passive 
linear optics and photodetectors with single photon sen- 
sitivity. If we apply the C-SHIFT gate to the qudits 1 
and k, where k = 2, . . . , M, then the GHZ state J23> is 
mapped onto product state of M photons 

\GHZ Md ) [X>^I 1 > 5 '1 ) I 1 ) 

d2---|l)dM- (28) 



D. State preparation with quantum memory 

The above described procedure for the preparation of 
the generic M-mode entangled state clearly shows that 
our approach is indeed general since it allows us, in prin- 
ciple, to generate an arbitrary state. But the success 
probability of the scheme will be exponentially small in 
general. The preparation procedure will therefore have 
to be repeated exponentially many times to produce the 
desired state. In the rest of this section we shall argue 
that the number of operations required to produce the 
desired state can be significantly decreased if a quantum 
memory is available. The situation is similar to that 
analyzed in Ref. 0] where the operations required to 
produce "cat states" was drastically decreased by using 
a recursive procedure that uses a quantum memory. 

The advantage of the quantum memory is that it allows 
one to wait until successive steps in the state generation 
protocol have succeeded before proceeding with the next 
steps. For instance in the scheme described in Fig. 5, 
one can first produce the states \fj) and store them in 
a memory before proceeding with the production of the 
entangled state (|19|l . Suppose that the state \f)j can 
be prepared with the probability Pj. Without quantum 
memory the probability of simultaneous preparation of 
the d states \f)j is P = Yi^i^i an d the number of 
required operations thus scales as NO oc 1/P. If, af- 
ter a successful preparation, we store each state |/) in a 
memory, we reduce the number of required operations to 

NO^U l l p i- 

A more elaborate use of a quantum memory is to 
build the state using a recursive procedure based on the 
Schmidt decomposition. The method is best explained on 
an explicit example. Suppose we would like to prepare 
a four-mode state \^)abcd ■ The Schmidt decomposi- 
tion of this four-mode state with respect to a bi-partite 
splitting into the AB and CD modes can be written as 



follows, 

d 

\iI>)abcd = 'J2\<Pj)ABWj)cD , (29) 

3=1 

where \tpj) ab and \itj)cD form orthogonal bases in some 
subspaces of the Hilbert space of modes AB and CD, 
respectively. The number of the terms in the Schmidt 
decomposition (|29|l is bounded by d < (N + l) 2 while, in 
general, (N+ 1) 4 terms were necessary in the general pro- 
cedure of Section IIIC. The decomposition (|2*9"|) suggests 
that we could obtain the state \iP)abcd via entanglement 
swapping if we first prepare two entangled states 

d 

\^)ABX =^2\(pj)AB\l)jX, 
3=1 
d 

\&)CDY = J2Mcv\l)iY- (30) 
3=1 

Now the entangled states \<&)abx and \$')cdy can be 
prepared using the general procedure of Sec. IIIC. How- 
ever, for states of the specific form H30|) . a simpler proce- 
dure can be devised. One first generates d states ab 
where the tilde indicates that \(pj) is related to \tpj) via 
relations similar to those described by Eqs. JUJ. The 
state \&)abx is then prepared with the help of the de- 
vice shown in Fig. 5 where the mode B of \<fij)AB is sent 
to the jth input port of the interferometer. The only 
modification is that we must also erase the spatial infor- 
mation carried by the A modes of the states \<pj) ab and 
perform the quantum state truncation after the erasing. 
Briefly, all d modes Aj are combined on an array of d — 1 
beam splitters and the detectors monitor the first d — 1 
output ports. If these detectors do not register any pho- 
ton, then the erasing procedure succeeded and we apply 
the quantum scissors to the state in the dth output port. 

Clearly, this recursive procedure can be extended to 
any number of modes. The resulting scheme resembles 
a tree structure involving repeated applications of the 
entangling scheme shown in Fig. 5 and the entangle- 
ment swapping that produces an 2m-mode state from two 
m-mode states entangled with auxiliary single photons. 
This scheme relying on quantum memory may possibly 
lead to a substantial improvement of the generation prob- 
ability, similarly to the case of the two-mode iV-photon 
cat states [l8j . 

IV. CONCLUSIONS 

In summary, the scheme proposed in the present paper 
is universal in that arbitrary multimode entangled states 
of light can be probabilistically generated. The resources 
required are passive linear optics, single photon sources, 
strong coherent states, and detectors with single-photon 
resolution. This result clearly illustrates the versatility 
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and power of the approach relying only on linear optics 
and single photons. 

It is fair to say though, that the suggested setup is 
rather complicated and involves several nontrivial opera- 
tions such as the quantum non-demolition measurement 
of a single photon, quantum scissors, and transformations 
of photonic qudits encoded as a state of a single photon 
in d spatial modes. We have shown that all these op- 
erations can be (probabilistically) implemented with the 
resources that we consider here. The resulting procedure 
may become very complex for general states. However, it 
should be emphasized that building a state like 12311 must 
necessarily be complicated because of the large number 
of parameters [(N+l) M — 1 complex numbers] that char- 
acterize this state and that must be fixed by the prepa- 
ration scheme. We have also argued that the use of a 
quantum memory can decrease the number of operations 
required to produce the desired state. 

We hope that this paper will stimulate the efforts to- 



wards experimental demonstrations of the basic building 
blocks of our scheme, such as the conditional generation 
of single-mode finite superpositions of Fock states and 
the preparation of two-mode entangled TV-photon states 
required for the quantum state truncation. 
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